3D crystal: how flat its flat facets are? 
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Abstract 

We investigate the hypothesis that the (random) crystal of the 
(— )-phase inside the (+)-phase of the 3D canonical Ising model has 
flat facets. We argue that it might need to be weakened, due to the 
possibility of formation of extra monolayer on a facet. We then prove 
this weaker hypothesis for the Solid-On-Solid model. 
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1 Introduction 



Consider the classical Ising model {as = ±1} , given by the Hamiltonian 

\s-t\=l 

in the cubic box Tat C Z'^ of side A^, with periodic boundary conditions and 
at the temperature which is low enough. Let us impose the canonical 
constraint: 

where |Tiv| is the volume of the box T^v, the constant 6 satisfies ^1 — {^Y^'^'^ < 

h < m , and where m {[3) is the spontaneous magnetization. Then the 
typical configuration a under the Gibbs canonical distribution will have a 
crystal: namely, it will have one large contour V [a) (which is a surface of 
codimension one, of linear size ~ and of volume ~ N'^), randomly lo- 
cated, such that inside F we will see the minus-phase, while outside V the 
plus-phase.^ 

We are interested in the geometry of the (random) crystal F (cr) . It is 
known that under a certain scaling the shape of F tends to a non-random 
limit. Namely, if one shifts F so that its center of mass will be at the origin, 
and then scales F by a factor of ^ in every direction, then in the limit N ^ oo 
the random surface -^F will approach the non-random surface (/3) , the 
well-known Wulff shape. (In fact, (P) depends also on 6, but since this 
dependence is just a linear scaling, we will omit it.) The meaning of the word 
"approach" depends on the dimension d. In dimension 2 the surface 14^2 (P) 
is just an analytic curve, and a question of convergence of F to W2 (/5) is 
studied in [DKS, IS] in great detail. Namely, with probability going to 1 as 
N ^ 00, one can shift the curve NW2 (P) in such a way that the contour F 
will lie inside the A'^^/^-ncighborhood of NW2 (/?) . In dimension 3 the known 
results hold in a weaker sense: one should pass from the configuration a 
to its integrated magnetization profile, which is a function M„ (x) G [—1,-1-1] 
on the unit torus T^. Then the distance between ^^^M^- (x) , properly 

^Of course we will have the separation of phases for all values of b. Our restriction 
excludes the case when the minus-phase arranges itself into a strip, wrapping around the 
torus. 
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shifted, and the signed characteristic function 21^1/3 (/?) — 1 of the inside of the 
surface goes to zero as ^ 00, see [CP, B, BIV]. Notice that on 

a suitable coarse grained scale, refined results on the stability of the Wulff 
crystal w.r.t. the Hausdorff distance were obtained in [BI]. 

Unlike the curve W2 (P) , the surface VF3 (P) is not analytic; moreover, it 
contains flat pieces - called facets - provided that the temperature is 
below certain critical temperature - called roughening temperature. It 
is known rigorously that > Tc (2) , see [BFL, BFM], where we denote by 
Tc {d) the critical temperature of the li-dimensional Ising model. It is an open 
question whether is equal to Tc (3) or is strictly less, as the common belief 
is. The shape of the facets of W3 is also given by the Wulff construction, 
see [M] or [SI], sect. 2.5. 

On the microscopic level, it was proven by Dobrushin [D] that at suffi- 
ciently low temperatures, rigid interfaces occur for some Gibbs measures with 
specific choices of boundary conditions. Therefore it is a natural question to 
ask, in which sense the flat facets observed in the macroscopic crystals and 
the microscopic rigid interfaces arc related. In this paper we want to discuss 
the question of whether or not the random crystals F {a) themselves have 
flat facets, for N large. Clearly, the results concerning the L^-convergence of 
r (cr) to W3 {j3) are perfectly consistent with either behavior. Some time ago 
one of us made the following conjecture, see [SI], sect. 3.4: 

Conjecture 1 (Probably wrong) Let the temperature j3~^ he low enough. 
Then the following event has probability approaching 1 as N ^ 00 : 
There exist six distinct 2D planes Lj = Lj (a) C T^v, i — 1, 6, two for each 
coordinate direction, such that the intersections Lj fl F (cr) are flat facets of 

F (a) . Namely, for every i 

i) diam (L^ n F {a)) > Ci {{3) diam (F {a)), with Ci {(3) ^^ifl as (3 ^ 00; 

[iTcL^inr^^^^^^^^ ^ ^2 m , with C2 iP) ^ 1/2 as P ^ 00, where by 
Area (Lj fl F (cr)) we mean the number of plaquettes o/F (cr) , belonging to the 
plane Lf, 

Hi) The asymptotic shape of the facets L^nF (cr) is given by the corresponding 
Wulff construction, see [M] or [SI], sect. 2.5. 

We beheve now that the above statement is a little bit too strong to be 
true. More precisely, it is almost true, except that one of the above 6 facets 
has an extra monolayer of the height one! So our reflned conjecture looks as 
follows: 
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Conjecture 2 (Hopefully correct) Let the temperature (5"^ he low enough. 

Then the following event has probability approaching 1 as N oo : 
There exist six distinct 2D planes Li = Li {a) C Tat, i = I, 6, two for each 
coordinate direction, such that the intersections fl F (cr) are flat facets of 
r (a) in the following sense: 



for every i diam {Li n T (cr)) > Ci (/?) diam (F {a)), with Ci {(5) ^^fljl 
as P —>■ oo; 



for every i except i — io — io (a) 
l^^XS- ^ ' ^''^ - 1/2 a. /3 - oo; 

Area(L,Q nr((T)) +Arca( (l^q +n,(, ) nr(o-)) 



. ^ — ■■ > C2 (P) , where is the unit vector 

[diam(LiQnr((7)jJ 

orthogonal to Li^ and pointing "away from" V (cr) . 

The meaning of the last statement is that on the facet Ljy fl V {a) there 
is another "monoatomic" layer of our crystal, having the shape (Lj^ + n^g) n 
r (cr) . The hmiting values ^^2/3 and 1/2 are coming from the fact that in 
the limit /3 — > 00 we expect F (cr) to approach the shape of the cube. 

At present we have no proof of this conjecture, and our paper is a result 
of an attempt to prove it. Namely, we prove here a weaker statement, and 
for a simpler - SOS - model. More precisely, we show that in the "canonical" 
SOS-model indeed a flat facet is formed, which may have an extra monolayer 
of particles. We formulate our result in the next section. In section 3 we 
further discuss it and we make various comments concerning the validity of 
the conjectures above. 

The heuristic explanation of our result is simple. Imagine that on one 
facet of the crystal F {a) we have two monolayers - the top one, Fi (cr) , 
located over the second one, F2 (cr) , with the size of the second one signifi- 
cantly smaller than the size of F (a) itself. Then we can enlarge F2 (cr) to the 
full size of the facet of F (a) , diminishing at the same time the monolayer 
Fl (cr) . It might even be that by that procedure the monolayer Fi (cr) will 
disappear completely. But in any case this procedure decreases the surface 
energy of the crystal. The reason for that is the same as for the fact that 
merging together two droplets into a larger one decreases the overall surface 
energy. In fact, we need here a slightly more general statement, that if the 
possible growth of the larger droplet is constrained by the container, then 
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still to grow it to the maximal possible size, while diminishing the smaller 
one correspondingly is energetically favorable. In the present paper it will be 
enough for us to have a zero-temperature analog of this statement, which is 
the content of the Lemma 8 below. This statement for the general case can 
be proven by the methods of the paper [SSI]. 

The second result of our paper deals with the question about the range 
of fluctuations of the random crystal around its limit shape. As we have 
said above, in the 3D case the known results about the closeness of the 
random crystal to its asymptotic shape are obtained only for the distance 
between them, while in the 2D case they are known to hold for the Hausdorff 
distance. Probably one cannot hope to extend this result to the 3D case 
at all subcritical temperatures. However it is reasonable to expect that such 
result does hold at very low temperatures. That was suggested already in the 
book [DKS] . Namely, though the solution of the Wulff variational problem 
is not stable in the Hausdorff distance, due to the possibility that thin long 
hairs can appear on the crystal, at low temperatures these hairs are highly 
improbable due to their energetic cost. Here we give an extra reason to 
believe it by proving "No Hairs" theorem, that for the low temperature SOS 
model in the N x N box the random surface fluctuates away from the flat 
facet by less than Cln A^, for some C < oo, and so the low temperature SOS 
crystal is always "clean-shaven" . 

We finish this introduction by pointing out the technical innovations of the 
present paper. Usually, to prove a result of such kind, one has to obtain the 
lower estimate on the probability of "nice" behavior of the random surface 
we are interested in, together with the upper estimate on the probability 
of its "ugly" behavior. Here the latter is easy, while the former is very 
hard, since this is the question about the typical behavior of the collection of 
contours which are strongly interacting, see [FPS]. We manage to establish 
our result by having only the upper estimate. This is both the strong and 
the weak point of our approach; we prove our theorem, but we do not have 
the technique to obtain the complete control over our model. 
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2 Statement of the Main Result 



Let (p = {ips G Z} be an integer valued random field, defined for s e Z^. Its 
distribution is defined by the "Solid-on-Solid" Hamiltonian 



\s-t\=l 



Namely, let A C be a finite box, |A| < oo, the configuration (boundary 
condition) tp. be given outside A, and the parameter /? > (inverse temper- 
ature) is fixed. Then we define the distribution Q/3,a,^ on the configurations 
O^K^ {^s ■ s e A} by 



Here 



<5/3,A,V- = exp { -I3H^ ((^ I V^) } / Z (/3, A, V^) . 
^a(<^|^)= l(<^VV'),-(<^VV')tl, 



(1) 



\s-t\=l 

{(f V i{j)t equals to (ft for t e A and to ipt for t ^ A, and the partition function 
Z (PjAjip) is a normalizing factor, making (1) a probability distribution. 

Our model of the crystal will be the distribution obtained from Q/3,a,^ by 
a suitable conditioning. Namely, we will consider the case when 

• A = An = {s : 1 < Si < N,i = 1,2} is a square N x N, 
. V'. = 0, 

• the volume constraint 

Viv(<^) = X<^. > AATS (2) 



is imposed, with A > fixed. 



We denote the conditional distribution (5/3,Ajv,V'=o (v 
P)3,N (95) • We do not keep A in this notation, since it will be fixed throughout 
the paper. We will use the notation Qj3^N for the unconditional distribution 
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<5/3,Ajv,V'=o (•) ■ We define the crystal C (</?) to be the body below the graph of 
V ■ 

C((^) = G : s G A,0 < /i <ip{s)] . 

To formulate our results about the facets we have to introduce the level 
sets. So for every (/? and every integer / > we denote by D {ip. I) the subset 
of all sites s in A, where (p (s) > I. We identify D {if, I) with the union of the 
closed unit squares centered at the corresponding points s. The connected 
components of the topological boundary of D (ip, I) will be called contours. 
The set of all contours will be denoted by A ((/?, /) . The sets D ((/?, /) can be 
disconnected; we denote by Di {ip, I) the collections of connected components 
of D{ip,l), i = 1,2,... which are mutually external. They will be called 
sections. By dDi ((/?, I) we denote the outer component of the boundary of 
the section Di((p,l). The set of contours dDi[ip,l) coincides with the set 
of external contours of the family A {(p>, I) . A section {if, I) will be called 
large, if 

\dDi{ip,l)\> KlnN, (3) 

where K is some big constant, to be chosen later. Otherwise it is called 
small. 

Consider now the level L — L{if) , which is defined to be the maximal 
value of l-s, satisfying the following condition: 

• \D {ip, l)\ > a {(3) N"^, where a (/?) is some small quantity, a (/?) — > as 
/5 ^ oo, to be defined later. 

Denote by Fi (ip) the level set D {ip,L {ip)) , and introduce also the no- 
tation Fj (ip) for the level sets D [if, L {if) — i + 1) . Our initial hypothesis 
was that the level set Fi {if) - the "First Facet" - is the facet sought, in the 
sense that |Fi ((/?)! > {1 — a {(3)) N'^ . However at present we cannot prove nor 
disprove this statement, and we think that it is not valid. In particular we 
cannot rule out the case of |Fi {if)\ ~ say. Still, we can show that a 

sharply localized jump of the function \D {ip, l)\ happens for typical ip-s: 

Theorem 3 Suppose the temperature is low enough. Then for the typical 
crystal the "Second Facet" is large: 

Pf3,N ■■ IF2 {ip)\>{l-a {P)) N^j^lasN^ 00, 

with some a {j3) —>■ as (3 ^ 00. 
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This result means that the crystal C (99) indeed has a horizontal facet, in 
the following sense: the level of height L ((/?) — 1 is almost filled with sites, 
whereas at the levels above L {(p) only few sites belong to the crystal. We 
do not know what happens at the level L {ip) , i.e. how big the First Facet 
really is. 

It is known that the SOS-model undergoes the roughening transition in 
temperature, see [FrSp]. At low temperatures the (unconstrained) SOS- 
model (without condition (2)) with zero boundary conditions is localized, 
while at high temperatures it diverges logarithmically with N. It is rea- 
sonable to conjecture that the roughening temperature Tf'-'^ is critical for 
our problem as well. In particular, it will mean that for every temperature 
> Tf^^, every £ > and for every pair m, n of integers 



for any value of the parameter a > 0, used in the definition of the level height 
L ((y?) . As our theorem shows, the behavior opposite to (4) takes place at low 
temperatures, and we conjecture that it is the case for all temperatures below 



To formulate the No Hairs theorem we introduce the boundary 9Fj (ip) 
of the i-th facet to be just the boundary dD {ip, L (ip) — i + 1). The theorem 
states that inside 9F2 (ip) the surface (p is almost flat, up to logarithmic 
excitations. 

Theorem 4 There exists such that for any (3 > Pq, one can find C > 
for which the following holds 

P^,jv{</? : 3s e Int (aF2 \ips - L{ip)\ > CluN} ^ Q as N ^ 00. 

3 Zero temperature Ising crystal 

In this section we discuss the relations between the Conjectures 1 and 2 
above. The Conjecture 2 is clearly a weaker statement, so it is not surprising 
that we can prove its SOS-counterpart, while we can not prove nor disprove 
the SOS-version of Conjecture 1. The real reason why the Conjecture 2 is 
simpler is the fact that it is valid at zero temperature, while the Conjecture 
1 is definitely not. 





r 
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The question about the shape of the crystal in the canonical Ising model 
becomes in the case of zero temperature the question about the isoperimetric 
problem in . Namely, we are looking into the following problem: let K be 
an integer, and we consider the family Vk of all subsets V <^1? containing 
precisely K sites. For every V we define the value \dV\ to be the number of 
plaquettes (of the dual lattice) in the boundary of V\ in other words, \dV\ is 
the area of the surface dV. We define Vk C Vk to be the subset consisting 
of minimal V-?,: 

V ^Vk^ \dy\ = min \dW\ . 

weVx 

In the following we will not distinguish the elements of Vk which differ by 
translation only, thus Vk becomes a finite set, so we can endow it with a 
uniform probability distribution. We want to take K to infinity and to look 
on the typical properties of the crystal shapes. However what we will see 
depends on the values of K. In case K — with M integer, the set Vk 
contains just one element, so the situation is trivial. To get some interesting 
behavior one has to choose the subsequence ^ oo in a special way. There 
are many different options here, and we describe just one of them. 

In the formulation of the theorem, which follows, the expressions "square 
with rounded corners" and "cube with rounded corners" are used. They 
mean here the following. Let k be an integer, and Yi, ...,^4 be four Young 
diagrams with the total number of cells less than k. Then the square k x k 
with four diagrams Yi,...,Y4 removed from its four corners is our "square 
with rounded corners". In the same way, a cube k x k x k with rounded 
corners is obtained from A;-cube by removing eight 3D Young diagrams (called 
also "skyscrapers") ^i, Ss, with the total volume below k, from its eight 
corners. We will call these diagrams as defects. 

Let < /X < 1 be a fixed number. We take 

K^^n^ + k{k-l) + l, 
where k = [iin] denote the integer part of iin. 

Theorem 5 As n ^ 00, a typical random shape V, drawn from the uniform 
distribution on Vk„, can be described as follows: 

V is a "cube with rounded corners" of size n, to one (random) face of 
which a monolayer is attached, which is a "square with rounded corners" of 
size k. These roundings have asymptotic shapes as n ^ 00: namely, letx>0 
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satisfies 



^2 2^13=^ (C (3)) 



2 



r6 



(so X is of the order ofk'^^^). Then each rounded corner of the square, scaled 
down by a factor x^^^, has asymptotic shape given by the Vershik curve, 

while each rounded comer of the cube, scaled down by a factor {k — 4a;) , 
has asymptotic shape given by the Cerf-Kenyon surface, see [CK], Theorems 
1.2 and 1.3. 

The curve (6) was obtained in [VKer]. The proof of the above result 
will appear later, see [S2]. The equation (5) is related with the asymptotic 
numbers of partitions and plane partitions of a large integer, see [SI], Sect. 
4.1 and 4.2. 



4 Proof of the Second Facet theorem 

We first prove a weaker statement, which, in fact, contains the main part of 
the proof of our result. Define £^(99) as the interior volume of all the external 
contours diD{ip, L{(p) — j + 1) of the level set D{ip, L{(p) — j + 1). 

Theorem 6 For any a > 0, there is such that 

V/3 > Po, P/3,N : E2 {^) >{l-a)N''}^lasN^ 00. 

Proof. The proof relies on energy estimates of the contours lying in the 
first and second facets. For a given height configuration ip, we denote by 
{'yi}i<Ki the set of all external contours of the family A ((/?, L{ip)) . These are 
just the external boundaries dDi{ip, L{ip)) of the connected components of 
Fi ((/?). Similarly, {Ti}i^K2 will refer to the external contours in F2 By 
construction the contours satisfy a compatibility condition, namely for any 
7i there exists Tj such that 7, lies inside Tj. 

We introduce two events; the first one, S, consists of configurations such 
that the volume contribution to Ei((/7) of external small contours of the first 
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2 ' 

tions for which the volume of the external large contours in the first facet is 

a 
2' 



facet is larger than ^N^, while the second, C, corresponds to the configura- 
lich the 

above |A^^, and also the volume of the external large contours in the second 



facet is smaller than (1 — a) A^^ : 
S 



L; |Int(7.)l>|iv4, (7) 

l> 7i small ) 

£ = L; Yl > E Il^t(r,)| < (1 - a)N' i . (8) 

^ 7i large Tj large j 

By construction, |Ei {(p)\ > |Fi {ip)\ > aN"^, so we can write 

Pp^N {if : |E2 {ip)\ < (1 - a) 7V2} < p^ ^^^S) + Pp^M . 

Thus, to complete the proof, it is enough to show that for /? large enough, 
there exists c = c(a, /?) > such that 

P0,n{C) < exp(-c7V) , (9) 
N ^ 



PpA^) < exp ( -c J ) . (10) 

The estimate (9) on the large contours will be obtained in Subsection 4.3 
and the estimate (10) on the small contours in Subsection 4.4. ■ 



4.1 A priory estimates on the height of the facet 

We start with very elementary estimates. Every SOS-surface ip e Qaj^ is 
made from N'^ horizontal plaquettes and a number of vertical ones; we de- 
note this last one by S {ip). Evidently, the distribution (3/3,Ar(<^) equals to 
exp{—pS {ip)} , up to normalization factor. Standard Peierls and counting 
arguments lead to the following simple estimate: 

Qp,m{^ ■.Siip) = S)< 3^'+^exp{-/35}. (11) 

Under Pp_N, the facet should be located with a high probability at a 
height of order A^. But for us a weaker statement will be sufficient. For K 
> 1 and A; > we define 

n = {(p; k< L{(p) < KN) . 
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Proposition 7 For any (5 large enough, any k fixed and K > K {(3) large 
enough 

Pp,N{n)>l-exp{-(3'N^), (12) 
where (3' diverges with (3. 

Proof. Our claim follows easily from (11) . Indeed, the property L{(p) > KN 

implies that 

S (if) > AK^/^N\ 

since for every level I below L{(p) we have | £)(</?, Z)| > a (/3) N"^. On the other 
hand, if the surface is such that all its sections ((/?, I) have the area below 
a ip) N'^ for / > then 

S ((^) > ^(^)^2 4\/H^iV = C {(3) N\ with C (/3) ^ oo as /3 ^ oo. 
Therefore 

since the denominator is always larger than exp(— /3AA^^). ■ 

4.2 Isoperimetric inequality - zero temperature case 

Here we prove the statement mentioned in the introduction, that merging two 
droplets together decreases the surface energy. More generally, just increasing 
the bigger one at the expense of the smaller one still makes the energy smaller. 
We prove here the corresponding statement for the 2D zero temperature Ising 
model only, while some generalizations are available by using the technique 
of [SSI]. 

For an integer V wc define L — L (V) to be the largest integer such 
that (y) < V, and we introduce r = r {V) = V — (V) . We denote by 
p = p (y ) the length of the shortest closed path on the lattice Z^, surrounding 
V unit plaquettes. Clearly, 

r 4L{V) ifV = L'^{V), 

p{V) = I AlIv) + 2 ifO<r{V)<L{V), 

[ 4L (V) +4 iiL{V)<r {V) < 2L {V) . 
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We will call p {V) the surface energy of the droplet V. In what follows we 
will identify the integers V with the collections of plaquettes from with 
perimeter p (V) , which will be called also droplets. Now, let Vi < V2 be two 
integers, and we suppose that for some and some (small) p > we have 

Vi > pN\ 

V2 < N'^. 

The second condition means that the larger droplet V2 fits inside the volume 
N X N, and the first one - that the smaller droplet is not too small. 
Let now D be any integer, satisfying the conditions 

Vi>D> pN\ 

Lemma 8 The transfer of the amount D from the droplet Vi to V2 decreases 
the total surface energy: there exists a constant x — x{p) > 0, such that 

(1 - X (p)) (p (Fi) + p {V2)) >p(Vi-D)+p (V2 + D) . 

Proof. We will show that the difference 

p{V^-D)+p {V2 + D)-p {V^) - p {V2) 

is of the order oipiVi) +p (V2) and negative. Since the function p [x) equals 
approximately to - more precisely, 

4v^ < P (a;) < 4v^ + 4, (13) 

- it is enough to show that the difference 

^/V^ - D + ^V2 + D - - ^2 
is of the order of \/Vi + and negative. Let us rewrite the difference as 




'VI 

and use the Taylor expansion of the function \/l + x. We get 




(14) 
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Now, since Vi < V2, the contents of all the odd brackets are negative, while 
the even coefficients are also negative, so the difference is negative as well. 
Finally, since all the values Vi, V2 and D are of the same order, the second 
term + y~/^ is of the order of ^/Vi + x/I^, and the proof follows. ■ 

4.3 Estimates on the large contours 

In this section wc will prove the estimate (9) . 

For a given integer £ > 2 and a compatible collection of large contours 
(7, r) = ({7i} , {^j}) we denote by </? ~ (7, F, the height configurations (p 
which satisfy: 

• the volume constraint Viv(<^) = "^seA^s ^ ATV^, 

• the only exterior large contours on the level sets Fi and F2 are given 

by (7,r). 

Then we have 



^ KN 

For a given triplet (7, F, i), we define the erasing map 

ip^(p= ((^^ = (^^_^ l{^gi„t(^.)} - ^ l{seint(rj)} J ■ (15) 
It maps injectively the set {ip ~ (7, F, i)} into the set 

I V^(^) > ATV^ - ^ |Int(7.)| - ■ 

Evidently, 

QpMv) = exp ( - /3 |7i| -PYI I^J'O QpM^) ■ 

i 3 
14 



Therefore we have the "Peierls estimate" 

KN 



Pp,M < Pp,^{W) -rY, E exp(-/?J]|7.|-/3^|r,|) (16) 

(.=2 (7,r)e/: i j 

Qp,n{ ^n{^) > AiV3 - |Int(7.)| - |Int(r,-)|) 



X 



The important quantity is the total volume of the interiors of the contours, 

thus we are going to average on all the possible contour shapes in order to 
retain only the information on the volume. Fix V G [|iV^, (1 — a)N'^] and 
consider the collection (7j) of the large contours such that 

Y.m^^)\^V. (17) 

i 

The optimal shape for a contour of volume V is a square of side length 
L — L{V) with (possibly) an additional layer of r (l^) sites such that 

V^L {Vf + r{V), r{V)e{0,..., 2L {V)} . 

In this case the following isoperimetric inequality holds uniformly over the 
collection (7^) which satisfy the volume constraint (17) : 

J2h\>^L{V). 

i 

Let 13' = 13 — 10. Summing over all the collections of contours such that 
Ei |Int(7i)| ^ V, we get 

J]exp ( -(3Y,\ii\) < exp(-/3'4L {V)) I E Il^^P ( " l^l^^l) 

{7i) i \(7i) i 

<exp(-4/?'L(V^)), 

where the final inequality is obtained by taking into account the entropy of 
a single large contour 



Y^\[e^v{-nii\)<\^+ E 3^exp(-10^)) 

(Ti) i \ i>K\nN J 



- 1 < 1 
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Plugging this inequality in (16), we get 

+ J2 Yl (2Li + 1)(2L2 + 1) exp ( - 4f3'{L, + L^)) 

y/^N<Li<L2<VT^N 

(18) 

Qp,n{ ^n{v) > XN^ - {Li + If - {L, + If) 



e=2 



X 



Qp,n{Yn('p)>XN^) 



Here we have indexed the volume Vi of the large contours in Fi((/?) by the 
parameter Li such that Vi = L\ + ri (with < ri < 2Li). Thus for a 
given Li, there is at most (2Li + 1) corresponding quasi-cubes. Similarly, 

the contours in F2 arc indexed by L2. 



The final step is to show that for \/^N < < L2 < \J\ — aN we have 

(19) 

This inequality combined with (18) will imply (9). 

For this we will use our Lemma 8, with D — min {(Li + if, N'^ — (L2 + 1)^} . 
From it we know that for /3 and N large enough 

4(1 - j) (Li + L2) > p ((Li + lf-D)+p ((L2 + If + D) . (20) 

As a consequence the function a (P) of the Theorem 3 must be chosen in such 
a way that 

X (a (/?))> 20//?. (21) 

Any choice of the function a {(3) as f3 ^ 00, consistent with (21) , is 
allowed in Theorem 3. So a {f3) vanishes as f] diverges, but nevertheless 
it cannot be too small. The heuristic reason is that for any finite /3 the 
macroscopic crystal has rounded edges; thus there exists a constant a' (/3) > 
such that it is no longer favorable to erase level sets which have the volume 
larger than (1 - a' {(3)) N^. 

Now to any height configuration (/? in {Vjv(v') ^ XN^ — {Li + if — (L2 + 
1)^}, we associate the configuration ip, defined by 

'Ps^ 'Ps + l{se[(L2+i)2+D]} + l{s6[(Li+i)2-D]} , Vs = (Si, S2) G A, (22) 
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where [n] denotes the square droplet n x n. Here one needs of course to fix 
the position of the smaller droplet [(Li + 1)^ — D] to be inside the larger one, 
[{L2 + 1)^ + D] . This however holds automatically, since either the droplet 
[(Li + If - D] is empty, or [(L2 + lf + D] = N x N. 

The correspondence (p (p maps injectively the set {Vat (99) > XN^ — 
(Li+1)2-(L2+1)2} into {Yn{<p) > XN^}- Furthermore, the energy difference 
between the height configurations (p and (p is bounded by p ((Li + 1)^ — D) + 
p((L2 + l)2 + D),so 

Qp,n{v : Viv(y^) > XN^ - {Li + 1)^ - (L2 + 1)^) 

< exp {P [p ((Li + iy-D)+p ((L2 + ly + D)]} . 
Combining this inequality with (20), we conclude that (19) holds. 



4.4 Estimates on the small contours 

This subsection contains the proof of the estimate (10) . 

We follow the scheme of the proof used to control the phase of small 
contours in the Ising model (see [SS2]). Define the subset A W = (i^^ x Z)^n 
A, where — 2K\nN. For any site s in A such that ||s||oo < -^^at, the shift 
of A'^°) by s is denoted by A^^\ To any collection of contours 7 = (7,4), the 
number of sites in A^^^ belonging to the interior of a contour in 7 is denoted 

If belongs to 5" and 7 is the collection of all the small contours in Fi(<^), 
then Yj\\s\\oo<Kn-^^ — i'^^ ■ Thus there exists at least one site s such that 

a / AT ^ ^ 



Denote by A* the set of collections of exterior small contours 7*^*^ — (7^*^), 
such that the condition (23) is fulfilled. Now for 7^*) e A* and i > 2 we 
introduce the sets S (7^*^-^) of the height configurations </?, which satisfy the 
properties: 

. L((^) = i, Y^i^) > XN^ 

• the contours from 7*^*^ are among the exterior contours of Fi(<^). 
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We get 

PMs)<pp,Nm + J2 E E E 

We proceed as in the previous subsection and erase all the small contours, 
belonging to the set 7'-*^. The total volume contribution of these contours is 
always smaller than N"^, so we get 

pp,N{s)<Pfi,Nm + T. E E e^p(-/5Ei^l'^i) (24) 

^=2 \\s\\oo<Kn 7(«)eA^ i 



X 



Using a shift of the height configurations by 1, we easily see that 

(3ftK(V„(43) > AJVS) - 1'^ ' 

In order to complete the derivation of (10), it is enough to prove that 
there is c > such that uniformly in s, 

J2 exp ( - 5^ |7f ) I) < exp (-ciVVn) ■ (25) 

The occurrence of a small contour surrounding a site io is bounded by 
^exp(-/3|r|) < 5^£3^exp(-/3£) = , 

where vanishes as P goes to infinity. 

Denote the number of sites N'^/Kjj in A^*) by Mjv. The small contours 
surrounding different sites of A*^*^ are independent, so we will obtain the 
bound (25) from an upper bound for a large deviation of a system of Mjv 
independent random variables. We have 

j: exp(-/3i:k'"i)< E 



fc>fMjv 
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where the last sum is the probabihty of the following event: 

Let ^1, ...,^Mn tie i.i.d. random variables, taking values 1 with probability 
q/3 and with probability 1 - q^. Then ^k>qMM {^k)'ipi^-9(3)'^'^~'' = 
Pr + ... + C^Mn — ^^n} ■ It is well known that such probability can be 
estimated from above by exp(— c^Mjv), where is a positive constant for 
/3 large enough. For a reference one can consult, for example. Lemma 10 
and Corollary 11 of [MRSV]. On the other hand — > as /3 — > oo, so the 
estimate (25) follows. 



4.5 End of the proof 

Thus far we have proven that the quantity E2(</?) - the area of the external 

contours {Ti}i^K2 of the second facet F2(v?) - is typically above the level 
(1 — a A^^. We are going to explain now that this in fact easily implies 
that the area of the facet F2{(p) itself has to be above the level (1 — 2a N^. 
Indeed, suppose that two events happen: 

E2 (if) > (1 - a (/?)) N^, and IF2 {ip)\ < (1 - 2a (/?)) N^. (26) 

We will show that its probability vanishes as A?" — > 00. 

To see this, let us introduce the second order external contours IVif , 

L J j<K2 

by defining them to be all the external contours of the collection 

A (99, L{(p) — 1) \ {^i}i<K2- (So the contours {rj}j<x2 should be called the 

first order external contours.) 
Under (26) we have that 

J2 |Int(f,)|>a7V2. 

j<K2 

Note also that the erasing map is now given by 

\ j<K2 

(compare with (15)), so 

> Yn{^) > AiVl (27) 
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As above, we split the second order external contours \ Tj> _ into 

I ' ij<K2 

small and large ones. We introduce two events; the first one, S, consists 
of configurations such that their small second order external contours of 
the second facet have a total volume larger than |iV^, while the second, C, 
corresponds to the configurations with the volume of the large second order 
external contours of the second facet is above ^N'^ : 



^; Yl |Int(f,)|>^iV2 

Tj large 



Then we can estimate the probabilities Pj3^i\f{C), Pg^N{S) by repeating the 
estimates for Pgj\f{C), P^^Ar(iS), obtained above. In fact, the corresponding 
estimates are even simpler, because the analogs of the estimates (18) , (24) , 

do not contain the factors '^^■^(^ Vjv(y)_A ) ^ 



5 Proof of the No Hairs theorem 

The excitations of the microscopic crystal around the second facet will be 
called hairs. There are two kinds of hairs: the up-hairs and the down-hairs. 
The up-hair of the SOS-surface </? is a sequence 71, . . . , 7if of contours, such 
that: 



^iedD{ip,L{ip) + i), (28) 

• contours 7^ are ordered by inclusion, which means that Int (7^+1) Q 
Int (7j) for alH > 1, 

• Int (71) C Int {dF2 {ip)) , 

• the sequence 71, . . . , 7ij is maximal, in the sense that there is no longer 
sequence of external contours, satisfying all of the above, of which our 
sequence is a subsequence. 
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We denote such an up-hair by F = {71, . . . , ^h}- 

The down-hair is defined as a similar sequence T = {71, ... ,7^} of con- 
tours, except that in (28) the sign is opposite: 7^ G dD {ip, L {(p) — i) . The 
value H will be called the length of the hair. 

Clearly, our statement is equivalent to proving that the probability of 
occurrence of a hair with length H > ClnN vanishes in the limit as — > 00. 
We denote such an event by Ac- 

In what follows we will treat only the up-hairs, since the case of the down- 
hairs is simpler, as there is no volume constraint (see the argument of the 
subsection 4.5); therefore, in the rest of this section we will use the term 
"hair" instead of "up-hair" . 

By our definitions we have that |Int(7i)| < aN'^. We introduce now the 
sequence 

Vr{N)^a—, r = 1,2,..., 

and we will characterize each hair by the amounts Hr{r) of its contours ji 
such that 

Vr+i{N) < |Int (7i) I < Vr{N) . (29) 

Naturally, we need only these Vr{N)-s which are > 1, so we define Rn to 
be the largest value of r such that the scale Vr{N) > 1. We fix also the 
intermediate scale R'j^ < Rn : 



Then for every r < i?^ we have 

Vr{N) > Ci{lnN)\ 

and the choice of the constant Ci is made in such a way that the contours of 
volume larger than Ci(ln A'^)^ are large contours, i.e. their perimeter is larger 
than KlnN. 

Finally we introduce the sequence hr{N), defined as follows: 



(4, li r < R!j 

hr{N) = I T-I'^N-^ In N, if r > R'j 



■N 1 
I 

N 1 



where C2 is chosen such that 

y/a C2 In ^/a 



hn'^ (N) = ' ^ =^C2>10. (30) 
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By definition, 

Rn Rn 
r=0 r=R'j, 

for some C3 > 0. 

Let us fix the constant C of our theorem to be much larger than C3. If 
the length of F exceeds C In N , then we define the value ro as the first index 
r for which the bound Hj.{T) < hr{N) is violated. Then define £ as the first 
index such that 

|Int(7^)| < tV„(iV). 

To summarize, to any F with length larger than ClnA^ we associate a pair 
(ro, i) and the subsequence of contours {7^, 7^+1, . . . , •je+hro-i} T, for which 
(29) holds. We denote by (ro, {7^, . . . , ^£+hr -1}) the class of all such hairs 
F. 

The strategy of the proof will be to apply a Peierls type estimate (under 
the volume constraint) to the section of the hair made of the exterior contours 
{7£, . . . ,'ye+hro-i}- 

Pp,n{Ac}< (31) 

Rn CalnN 

< ^ ^ ^ PfS^N {(p has a hair F in the class (ro, i, {7^, • . . , 'yt+hr^-i}) } , 

ro=l 1=1 {7i} 

where the sum is over the collections of contours {7^, . . . , ^i+K^^-i}- In order 
to estimate (31), two cases have to be distinguished. Either tq is smaller 
than R'j^ and all the contours 7^, . . . , 7^_|_/j^^_i are large, or ro is larger than 
i?'^ and one has to rely on more delicate estimates, taking into account the 
fact that these contours might be small. 

Case 1 : ro < -R^- 

In this case /i^o — 4 and we have 

3 

^|Int(7,+OI <4^;.oW, 

i=0 
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J2 > 4(4V^;.o+l(^)) • 



1=0 

Applying the Peierls estimate as in the proof of Theorem 6 (see relation 
(15)), we get for every ro — [0, R'j^] and every height ^ < Csln 

^P^,^r{re(ro,£,{7^,...,W)} (32) 

in} 

' ^ ^iiv(Vjv(v^)>AAr3-4t;,„(7V)) 



< ^exp I -/3j2he+i\ 

{li} \ i=0 



Qp,N (Viv(<^) > AiV3) 



The loss of volume can be compensated by adding a single square contour 
of side length A[y^ Avro{N)] + 1, in the same manner as in relation (22) above. 
Furthermore the entropy of the four large contours can be easily bounded, 
so that we obtain 

R'n CslnAT 

E E E^A^{r^(ro,£,{7^,...,7m})} (33) 

ro=0 £=Q {7,} 

R'n C3I11JV 

<E E (-16(/9 - ln3)^yVr,+i{N) + 4p^4vr,{N)^ 

ro=0 e=o 
< exp{-Ci3lnN) , 

where Cp is a positive constant for (3 large enough. 
Case 2 : To > R'j^. 

First notice that for any collection of contours in the r.h.s. of (31), we 
have for the area 

i+hrQ-l 

^ |Int(7,)| <^;.o(^)^ro(iV), 

i=e 

while for the boundary we get 



E hi\>Wvro+l{N)hr,{N) 
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Since 

J\f2 (~i 



for r > i?^ = In (^c0^) we have 



vJN)hJN) < , ^ NlnN = J^MlaNf 



Ci(lnJV)2 

On the other hand, 



V^W^^r-(A^) = y^^C-s 2^/'N-' In AT = ^1^2 In TV. (34) 
Thus the energy Yfi^e'^° ^ the collection of contours exceeds 

and therefore, due to (30) , is much larger than the one of a single square 
contour with the same area, since its energy equals 



A^yVr,{N)hr,{N) < AsJ V^C2 In N. 

Nevertheless these contours can be small, and one has to estimate their en- 
tropy carefully. 

The number A/'(L£, . . . , L^+h^^_i) of compatible contours {7^, . . . , 7£+/,^^_i} 
with respective length {L^, . . . , can be estimated by 

J\f{Le, W„_i) < 3^^ H {vr,{N) 3^') . 

i=e+i 

To see this, one chooses first the contour 7^ of length Li at a random position 
in the box, then the other contours pile up above it, so their entropy is given 
by the second term, where Vro{N) is the maximal area of 7^. 
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We proceed as in (32) and get for a given tq and height £ 



{7i} 



< ^cxp 

{7i} 



V 



ne+i\ 



i=0 



Ni'VNi^)>^N^-Vro{N)hr,{N)) 



<5^iV2^xp K,{N)\n{vr,{N)) + {\n3-P) ^ + 4/?^AU^V)M^ | , 

{Li} \ i=i 



where each Li ranges in Vro+i{N) , Vro{N)], and so 



Summing over the sequence {L^}, we see that 
J2 Pl3,N {r ~ (ro, ^, {7£, . . . , -ft+hr^-i}) } 



(35) 



{7i} 



< iV^exp (Jgg - - 10) + 

Recall that K^iN) > hR'^{N) > 10 (see (30)), /iro(^)\/?V„+i(iV) 
yfC.toiV (^e (34)), while 
indices Tq and we obtain 



< 1. Thus summing (35) over the 



Rn Cain at 

E E E^/^.^{r~(^o,^,{7^,...,7m.o-i})} 

-Rjv Cain AT 

< 

< exp(-C^ln7V) 



E;E"-^exp(;fc.ln;v(l-4(,-10, + ^)) 



where C/j is positive constant for /? large enough. 
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Combining this with (33), we conclude the proof. 
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